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Abstract. We present the first systematic spectro-polarimetric study of Lu- 
Zj ' minous Blue Variables (LBVs), and find that at least half those objects studied 

, display evidence for intrinsic polarization - a signature of significant inhomo- 

geneity at the base of the wind. Furthermore, multi-epoch observations reveal 
t ' that the polarization is variable in both strength and position angle. This ev- 

idence points away from a simple axi-symmetric wind structure a la the B[e] 
supergiants, and instead suggests a wind consisting of localised density enhance- 
ments, or 'clumps'. We show with an analytical model that, in order to produce 
the observed variability, the clumps must be large, produced at or below the pho- 
tosphere, and ejected on timescales of days. More details of LBV wind-clumping 
will be determined through further analysis of the model and a polarimetric 
. monitoring campaign. 
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1. Introduction 



Luminous Blue Variables (LBVs, iHumphrevs fc Davidsoi3ll994L & Nota these 
' proceedings) are very massive evolved stars that are found right at the top of 

the HR diagram, roughly the same region as the B[e] supergiants. Like the 
B[e]SGs, their spectra contain many emission lines due to their strong winds. 
Unlike the B[e]SGs, they are strongly photometrically variable, with amplitudes 
of 1-2 mags over timescales of a few months to years. The amount of material 
ejected by the LBVs is huge; their mass-loss rates are among the highest known 
(KT 5 - 10~ 4 M^yr- 1 ) whilst they c an th row off many solar masses during giant 
eruptions (e.g. n Car, iMorris et alJll99 Sh. The LBV phase therefore represents 
a crucial stage in a massive star's evolution. 

The products of the extreme mass-losin g episodes of LBV s can be seen in 
the form of their surrounding nebulae (e.g. iNota et alJll995f h However, it is 
unclear whether the bi-polar nature of these nebulae is due to a pre-existing 
density contrast, or if the wind itself is axi-symmetric. Such a wind may be 
linked to the asphericity of supernova explosions and the beaming of gamma- 
ray bursts; whilst wind asphericity is of major importance to stellar evolution , 
as it can lead to large over-estimates of mass-loss rates (e.g. iBouret et al.l l2005). 

The study of the inner-wind morphology of these stars is not straight- 
forward. They are too far away to directly image the inner wind, and spec- 
troscopy alone gives no unambiguous geometry information. The only tool ca- 
pable of determining the present-day mass-loss geometry is spectropolarimetry, 
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and we present here a comprehensive spectropolarimetric study of all LBVs in 
the Galaxy and the Magellanic Clouds. 

1.1. Introduction to Spectropolarimetry 

The technique is based on the following: the continuum emission of a star is 
electron-scattered by the ionised material at the base of the wind, while the 
line-emission, which forms in the wind over a much larger volume, remains 
essentially unscattered. If the material at the base of the wind is aspherical 
(e.g. an equatorially-enhanced flow), the continuum light will be left with a net 
polarization perpendicular to the plane of scattering. The line-emission, which 
is unscattered, remains unpolarized. The dilution of the polarized flux by the 
line-emission can be seen as a drop in polarization across the emission line (see 
Fig. m>. As the majority of the polarization occurs within a couple of stellar 
radii (|Cassineili et al.lll987n . by looking for changes in polarization across strong 
emission lines such as Ha we can look for inhomogeneous structure at the base 
of the wind. 
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Figure 1. Polarization spectra of two LBVs, AG Car (left), and HR Car 
(right). The bottom panels show the intensity spectrum in the region of Ha; 
the middle panel and upper panels show the degree of polarization and PA 
as a function of wavelength respectively. The changes in polarization across 
the stars' e mission lines are indi cative of aspherical structure low in the wind. 
Taken from lDavies et alJ l|200ifl . 



2. Spectropolarimetric observations of LBVs 

Using the above-described technique, we have observed all known LBVs in the 
Galaxy and Magellanic Clouds. We find that at least half of those observed 
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show changes in polarization across Ha, and therefore signatures of asphericity 
at the base of the wind. We put this 50% detection rate as a lower limit, due to 
the difficulty in achieving the target S/N of 1000 (0.1% precision) for the fainter 
stars. This compares with detection rates of 20% and 25% in similar studies of 
Wolf-Rayet stars and O supergiants respective ly, where positive detectio ns were 
very much the exception rather than the rule (|Harries et alJll998t I2002T ). 



p 

o 

li' 

O 

i ^Interstellar 

o \ 

o 

0.2 0.4 0.6 0.8 

% Q 0.0 0.5 1 .0 

%Q 

Figure 2. Left: Polarization vector diagram of P Cyg. The magnitude of 
the vectors between the continuum measurements (open circles) and the zero- 
point polarization of P Cyg (marked 'interstellar') represent the strength of 
the intrinsic polarization; while the angle of the vector with the Q axis rep- 
resents the polarization's position angle (PA). The polarization is variable in 
strength and PA on timescales of days. Taken from INordsicc k et al.l (|2001). 

Right: Same as the figure on the left but for AG Car. The zero-point polariza- 
tion (marked 'interstellar') h as been determined f rom the Ha depolarizations 
(dotted lines). Adapted from lDavies et all $2005). 

In addition, we find that intrinsic polarization is more likely to be found in 
those objects with strong line-e mission, and those t hat have been variable by ^1 
mag in the last ~10 years (see IPavies et al~ll2005h . This may be linked to the 
fact that the region occupied in the HR-diagram by the LBVs is close to certain 
instabilities, namely the bi-stability jumps and the Humphreys-Davidson limit. 
Strong variability in this regime means that they can often appear to stray into 
these zones, possibly resulting in erratic mass-loss behaviour, leading to strong 
line-emission and wind-inhomogeneity. 

2.1. Temporal Variability 

For four of our stars (AG Car, HR Car, P Cyg and R127) there exists multi- 
epoch data in the li t erature and archives dSchul tc-Ladbcck et alJ Il994l. 1 1993; 
Leitherer et alJ Il994l : IClampin et alJ Il995l : iTavlor et alJ Il99ll : iNordsieck et alJ 
200lh . For all four stars, the intrinsic polarization they exhibit is variable in both 
strength and PA (see Fig. |2J). This is inconsistent with a steady axi-symmetric 
wind, as observed in Be stars, where we would expect the PA to remain constant. 
To explain the variable polarization, we considered four possible scenarios, which 
are summarized belo w (see Fig. Efl. Fo r a full discussion on each explanation 
and their merits, see lDavies et al.1 (|2005h . 
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• An axi-symmetric wind with a variable optical depth, causing the polarized 
emission to originate from different latitudes as a function of time. In 
this scenario, polarization at any epoch would be aligned with one of two 
perpendicular directions (see Fig. As this is not observed, this rules 
out this scenario as the sole cause of the variability. 

• A 'flip-flopping' wind, where the wind's enhancement-plane switches be- 
tween equatorial and bi-polar. This would again produce polarization in 
two perpendicular planes and can be ruled out. 

• Binarity, similar to the mecha nism producing t he po larization observed in 
WR-0 binary systems, see e.g. lSt.-Louis et all ( 1987). This is unlikely due 



to the physical properties required of the companion; it would need to have 
comparable luminosity and/or mass-loss rate to the LBV. It is unlikely 
such an object would have not already been detected through doppler- 
wobble, composite spectrum, or X-rays from the wind-wind collision zone. 

• A 'clumpy' wind, where the polarization is caused by light scattering off 
localised density-enhancements low in the wind. As the clumps move away 
from the star in the wind and new clumps are constantly being formed, 
the distribution of clumps and hence the net polarization will appear very 
different from epoch to epoch. This would produce the stochastic polar- 
ization variations we observe in LBVs (cf. Fig. I^J), and is the most likely 
explanation. 

A clumpy in ner wind is consistent with the clumpy n ature of these stars' in- 
ner nebulae (e.g. iNota et ai", Hl99fit lCh esneau et al.ll2000l) and the spectro scopic 
variability of their discrete absorption components (e.s. IStahl et al.ll200ll ^. It is 
worth noting that wind-clumping is also the favoured explan ation of similar po- 
larimetric behaviour seen in WRs by e.g. lRobert et all (|l989l ). who also find that 
the amplitude of the variability decreases with increasing terminal wind velocity. 
This is consistent with the LBVs, which have much slower winds (~200kms _1 , 
compared with ~1500kms~ 1 ), with greater polarimetric amplitude (~1% com- 
pared with £0.2%). 

As an aside, the first few polarimetric observations of AG Car bv lSchulte-Ladbeck et all 

1 1994 1 suggested some form of axi-symmetry, leading the authors to suggest one 



of the first two explan ations. However, these scenarios are not s upported by sub- 
sequent observations ( Leitherer et al.l|l994l : iDavies et al.ll2005t ). This highlights 



the pit-falls in drawing conclusions from single-epoch polarimetry data. 



3. Simulations of clumpy winds 

In order to begin to get a handle on typical clump parameters, we have con- 
structed an analytical model which simulates the polarimetric variability due 
to wind clumping. Clumps are given a typical size, optical depth, and ejection 
timescale. Once ejected radially from a random position on the surface of the 
star, the polarization of one clump is tracked as it accelarates with a /3 = 1 
velocity law. The polarization due to all the clumps in the wind is calculated as 
a function of time, as new clumps are ejected and move through the wind. A 
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a) 



b) 



c) 




Figure 3. Cartoons of three different scenarios capable of producing variable 
polarization (left), and the corresponding Q — U behaviour we would expect 
(right). Symbols are the same as in Fig. 

(a) : an axi-symmetric wind. When the PA of polarization is perpendicular 
to the disk, we measure polarization along one vector in Q — U space. If the 
PA were to rotate by 90° due to e.g. an opacity increase at the equator or a 
flip from equatorial to bi-polar wind, we would measure polarization along a 
vector on the opposite quadrant of Q — U space. 

(b) : binarity, e.g. a WR-0 system. The O star's light is polarized due to 
scattering off the wind- wind collision zone. As the density enhancement region 
rotates with the binary system, the polarization describes a double-loop in 
Q — U space. 

(c) : wind-clumping. Each clump in the wind scatters the starlight. If the 
clump distribution is not spherically symmetric, a net polarization results. As 
the clump distribution changes with time, so does the polarization. 
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Figure 4. Simulations of the time-dependent polarization produced by a 
random distribution of clumps in a spherical wind. The simulations use typical 
LBV parameters of M — 30M Q , B = 3OOR.0, and Voo =200 km s , and 
the clumps have r e = 1 and ejection timescale of 7 days. Runs using three 
different clump sizes of 0.4R*, 0.3R*, and 0.2R* are shown. From Davies 
(2005, in prep). 



typical run is shown in Fig. 01 Full analysis of the models is pending (Davies et 
al. 2005 in prep), but initial results are summarized below. 

It can be seen that, in order to produce the observed levels of polarization 
and polarimetric variability, the clumps must be large (^ 0.2R+) and ejected 
every few days. Increasing the density of the clumps (and hence the number of 
scatterers) does not allow the clump size to be reduced, as Monte-Carlo stud- 
ies have shown that the pola rization per clump plateaus at the r e ^ 1 level 
(|Rodrigues Magalhaesl l2000D . Also, the polarization per clump falls off dra- 
matically with distance. This means that (a) the clumps must be produced at 
or below the photosphere, no t condense out of the wind at larger distances due 
to radiative instabilities (e.g. iRunacres fc OwockilEiol ^l: and (b) there must be 
a constant supply of clumps. 

For a gi ven mass-loss rate and density contrast between clump and ambient 
wind (~ 20, iNordsieck et al.ll200ll ). the polarization as a function of clump size 
and ejection timescale can be simulated. Through polarimetric monitoring of 
the LBVs we will constrain the clump ejection timescale through the signature 
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of jumps in polarization that ejection events cause (Fig. 0]). From this, we can 
ultimately determine the volume fil ling-factor of the c lumps, the parameter to 



which wind models are so sensitive ((Bouret et al.l 2005) 



4. Summary 

In the first comprehensive spectropolarimetric study of LBVs, we find that at 
least 50% show evidence of intrinsic polarization - over double the detection 
rate of WRs and O supergiants. Moreover, the polarization is stochastically 
variable, suggesting that the polarization is not due to axi-symmetric wind struc- 
ture but instead caused by significant wind-clumping. Preliminary modelling of 
this shows that the clumps must be very large, and produced at/below the 
photosphere - clumps condensing out of the wind at large distances cannot pro- 
duce the observed variability. Polarimetric monitoring will allow us to constrain 
further typical clump parameters, and ultimately determine the wind clump 
filling- factor, enabling more accurate mass-loss rate estimates. 
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